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The nature of the viscosity operative in hot, two-temperature accretion disks around AGN has been a long- 
standing, unsolved problem. It has been previously suggested that protons, in conjunction with the turbulent 
magnetic field that is likely to exist in the accretion disk, might be crucial in providing this viscosity. Several 
authors have recently determined diffusion coefficients for charged particles (cosmic rays) propagating in tur- 
bulent magnetic fields by means of extensive Monte Carlo simulations. We use the diffusion coefficients for 
protons determined by these simulations to find the effective mean free path for protons in hot accretion disks. 
This in turn yields good estimates of the viscosity due to energetic protons embedded in the turbulent mag- 
netic field of a hot, two-temperature accretion disk. We combine this with a simple two-temperature accretion 
disk model to determine the Shakura-Sunyaev a viscosity parameter arising out of this mechanism. We find 
that protons diffusing in the turbulent magentic field embedded in a hot accretion disk provide a physically 
reasonable source of viscosity in hot accretion disks around AGN. 

1. Introduction 

The microphysical viscosity mechanism operative in hot, two-temperature accretion disks around AGN has 
been a subject of intense research, and is a yet-unsolved question. In the steady state, such a viscosity mech- 
anism enables matter around the central black hole to lose its angular momentum and accrete onto it. Recent 
advances in the theory of the magnetorotational instability in cold accretion disks (Balbus & Hawley 1998 
and references therein) suggest the existence of a chaotic, tangled magnetic field, together with a relatively 
well-ordered large-scale toroidal field embedded in the accretion disk around a black hole. The "magnetic 
viscosity" arising from stresses related with the tangled magnetic fields is thought to be the operative viscosity 
mechanism in cold accretion disks. However, in hot, two-temperature accretion disks, where the protons are 
hot (Ti ~ 10 10 -10 12 K) and their mean-free path is comparable to (or larger than) the scale height of the disk, 
this picture is inadequate. The theory of MRI-like instabilities under conditions similar to those prevalent in 
hot accretion disks is still in its infancy (Balbus 2004; Quataert et al. 2002; Sharma et al. 2003). However, 
these initial developments do suggest that such an instability will also culminate in a tangled magnetic field 
embedded in the accretion disk. Assuming the presence of such a tangled magnetic field in a hot accretion disk, 
Subramanian, Becker & Kafatos (1996) (SBK96 from now on) showed that the operative viscosity mechanism 
was a "hybrid" one, arising from hot protons colliding with kinks in magnetic field lines. SBK96's treatment 
had £ = X co h/H, as a free parameter, where A co h is the coherence length of the magnetic field and H is the 
disk scale height. The quantity £ thus provides a measure of the extent to which the magnetic field is tangled, in 
comparison with the disk height. SBK96 were only able to treat one value of £ at a time, and could not consider 
the overall effect arising from the presence of a spectrum of £s (i.e., a spectrum of scales over which the field 
is tangled). Furthermore, they did not consider the presence of an ordered, large-scale toroidal field, which 
is suggested by MRI simulations, albeit in cold accretion disks. In the present treatment, we use results from 
extensive Monte-Carlo simulations of cosmic rays propagating through tangled magnetic fields published in 
the literature in order to remedy these shortcomings and arrive at a more rigorous estimate of hybrid viscosity. 
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2. Cosmic ray transport through turbulent magnetic fields 

Several authors have studied the diffusion and drift of charged particles across turbulent magnetic fields. While 
analytical perturbative studies in the low-turbulence regime have been carried out for a long time (e.g., Parker 
1965; Jokipii 1966) numerical experiments suitable to the high turbulence regime are relatively recent (Gi- 
acalone & Jokipii 1999; Casse et al 2002; Candia & Roulet 2004). There is a large scale magnetic field 
typically assumed to be present, in addition to the turbulent conponent, and such studies strive to characterize 
the motion of cosmic rays through turbulent magnetic fields by means of diffusion coefficients perpendicular 
and parallel to this large scale field. We follow the treatment of Candia & Roulet (2004), who have performed 
Monte Carlo simulations following cosmic ray trajectories in the presence of a tangled magnetic field for a 
variety of turbulence levels and magnetic rigidities. They provide convenient analytical fits to the parallel and 
perpendicular diffusion coefficients arising from their simulations. 



3. Hot protons diffusing through tangled magnetic fields: hybrid viscosity 

We use the analytical fits to the diffusion coefficients provided by Candia & Roulet (2004) to compute the 
viscosity arising out of hot protons diffusing through tangled magnetic fields. The advantage of this approach 
is that it takes the interaction of the hot protons with all the spatial scales of turbulent magnetic fields, and not 
only with one scale, as in SBK96. As suggested by simulations of the magnetorotational instability in cold 
accretion disks, we presume that there exists a large-scale magnetic field in the toroidal direction embedded 
in the accretion disk, in addition to a tangled field. We are interested in diffusion coefficients (and consequent 
viscous angular momentum transport) in the radial direction, which is perpendicular to this large-scale toroidal 
magnetic field. The perpendicular diffusion coefficient is, by definition, the one which characterizes motion 
perpendicular to the large scale, ordered magnetic field. We therefore use expressions for the perpendicular 
diffusion coefficient D± provided by Candia & Roulet (2004). 

From Eqs (18) and (19) of Candia & Roulet (2004), adopting the low rigidity assumption (p < 0.2) for 
simplicity, we get 

where we have used the disk height H for the maximum scale length L max - The quantity u rms = y/ZkBTi/rrii 
(where ks is Boltzmann's constant, Tj is the proton temperature and mi denotes the proton mass) is the proton 
rms velocity. We use u rms instead of the speed of light in defining D± as Candia & Roulet (2004) do, since the 
hot (Ti ~ 10 10 K) protons we consider have energies of the order of 1 Mev, and are therefore non-relativistic. 
Candia & Roulet's (2004) results are quite valid for protons in this energy range too, unless the proton thermal 
velocities are close to the Alfven speed, in which case coupling to Alfven modes needs to be taken into account. 
Although simulations of the nonlinear stage of the MRI in hot accretion disks have not yet been performed, 
corresponding simulations for conditions corresponding to cold disks reveal that the final plasma (5 saturates at 
values ~ 500. Under these conditions, the particle thermal velocities are far in excess of the Alfven velocity. 
The quantity p denotes the rigidity, which is defined as the ratio of the Larmor radius to the maximum scale 
length. The quantity a 2 denotes the ratio of the energy density in the turbulent magnetic field to that in the 
large scale magnetic field. The values of the parameters a±, N\\, , N± are given in Table 1 of Candia & Roulet 
(2004) for various kinds of turbulence. For the sake of concreteness, we use values corresponding to the first 
line of that table; i.e., Kraichnan turbulence. 
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The coefficient of dynamic viscosity is usually defined as i] (g cm -1 s _1 ) = NmvX, where N and m are the 
density and mass of the relevant particles, v is the relevant velocity (usually the thermal rms velocity) and 
A is the relevant mean free path (Spitzer 1962; Mihalas & Mihalas 1984). The diffusion co-efficient D is 
usually defined as D (cm 2 s" 1 ) = vX, where v is the relevant velocity and A is the relevant mean free path. 
Combining these two expressions, we get 77 = NmD. The coefficient of dynamic viscosity arising out of the 
radial transport of hot protons (i.e., perpendicular to the large scale toridal field) can therefore be written as 

r] hyh = N t m p D ± = 7.54 x 10 10 L> c r e;5 grants -1 , (2) 

where we have used Eq ([0 and Eq. (A6) of SBK96. The quantity Ni is the proton number density and m,i is 
the mass of a proton. 

Using Eq and Eqs. (Al), (A4), (A6), (A9) of SBK96 in Eq. (3.4) of the same paper, we get 

1 + r es = 3.07 x 10- 3 (J^) /r 1/2 f^yRl' 2 D c (3) 

and 

a hyb = 2.557 /r 1/2 £> c . (4) 

We have used the same notations as in SBK96; r es denotes the electron scattering optical depth, ahyb is the 
Shakura-Sunyaev viscosity parameter (Shakura & Sunyaev 1973) due to protons diffusing through tangled 
magnetic fields, M* is the accretion rate in units of solar masses per year, Ms is the black hole mass in units of 
10 8 solar masses and fx, /a and f^ are the relativistic correction factors used in SBK96 and references therein. 
The results for ahyb an d T es from Eqs © and © can be used in Eqs. (A7)-(A9) of SBK96 to determine the 
disk structure as a function of radius for a given accretion rate For a given set of parameters characterizing 
the turbulence (which yields a given value of D c via Eq.[T), Eq (fi} shows that the requirement r es > implies 
an upper bound on the accretion rate . 

In Figure 1, we have shown the results for the following set of parameters characterizing the turbulence: 
7 = 3/2, JVjj = 2, p\\ = 0.22, N± = 0.019, a± = 1.37. These values are taken from the first line of Table 1 of 
Candia & Roulet 2004, and characeterize a Kraichnan turbulent spectrum. The value of the magnetic rigidity 
p is taken to be 0.2 and the parameter characterizing the level of turbulence, a 2 , is taken to be 0.1. In other 
words, the proton larmor radius is 0.2 times the disk height, and the energy in the turbulent magnetic fields 
is 0.1 times that in the ordered magnetic field. The accretion rate = 1-6 x 10~ 4 , which is close to the 

D Ms 

allowed upper limit above which the optical depth tau es will become negative. This value of ^f- corresponds 
to M- = 7.27 x 10~ 4 , where Afg is the Eddington accretion rate. The solid line shows the value of «hyb> 
while the dotted line shows the quantity 0.1 H/R. The dashed line shows the proton temperature in units of 
10 12 K. We have found that the values of «h y b are in the range of 0.01 0.02 for a 2 = 0.5. 

4. Conclusions 

We have used results from extensive Monte Carlo simulations of cosmic ray transport through turbulent mag- 
netic fields to obtain a rigorous estimate for the hybrid viscosity in hot, two-temperature accretion disks. We 
find that the value of the Shakura-Sunyaev viscosity parameter arising from this mechanism is in the range 
0.01 > ahyb 0.02 for weak to moderate Kraichnan turbulence (0.1 > a 2 > 0.5). We have thus demon- 
strated that hybrid viscosity is a viable agent of steady state angular momentum transport in hot accretion 
disks. In SBK96, we compared the values of ah y b with those arising from pure magnetic viscosity a mag . In 
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Figure 1. Results corresponding to an accretion rate of M/Me = 7.27 x 10 4 . The solid line is a plot of the a parameter, 
while the dotted line shows the quantity 0.1 H/R. The dashed line shows the the proton temperature Tj in units of 10 12 K. 

that paper, we had restricted ourselves to a considering a single value of £ = \ co h/H at a time, and found that 
"mag = £ 2 - This facilitated a ready comparison between the two kinds of as. Here, however, we do not restrict 
ourselves to a single value of £, and use the entire turbulent spectrum to calculate ahyb- We cannot, therefore, 
resort to such a direct comparison. Simulations of the magnetorotational instability under conditions appro- 
priate for cold accretion disks suggest that a mag ~ 0.01, but it is not clear what role pure magnetic viscosity 
will play in hot accretion disks. The answer will depend crucially on the value of the plasma (3 attained in the 
steady state. 
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